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Abstract 

We provide detailed estimates of the proton lifetime in the context of simple su- 
persymmetric SU(5) grand unified models with an extra compact spatial dimen- 
sion, described by the orbifold /{Z2 x Zg) and by a large compactification scale 
Mc ~ 10^^ 10^^ GeV. We focus on a class of models where the grand unified 
symmetry is broken by the compactification mechanism and where baryon violation 
proceeds mainly through gauge vector boson exchange so that the proton lifetime 
scales as M^. We carefully compute Mc from a next-to-leading analysis of gauge 
coupling unification and we find that Mc can only be predicted up to an overall 
factor 10^^. The simplest model, where the dominant decay mode is vr'^e"^ and has 
no flavour suppression, is strongly constrained by existing data, but not totally ruled 
out. We also analyze models where some of the matter fields are localized in the 
extra space and proton decay is flavour suppressed. In models associated to anarchy 
in the neutrino sector the preferred decay channel is K^v and the lifetime can be 
within the reach of the next generation of experiments. 
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1 Introduction 



Grand unification can be considered a cornerstone in our quest for unification of particle 
interactions. Many properties of tlie standard model of electroweak and strong interac- 
tions that seem mysterious or accidental, like the particle content, the cancellation of gauge 
anomalies, the quantization of the electric charge, appear natural in the context of grand 
unified theories (GUTs). The quantitative success of gauge coupling unification in super- 
symmetric (SUSY) GUTs is remarkable and the vicinity of the GUT scale to the Planck 
scale is quite intriguing. Also the observed smallness of neutrino masses, suggesting the 
breaking of the total lepton number L at very high energies, supports the idea of a new 
threshold in particle physics at the GUT scale. 

Despite the attractiveness of the GUT idea and its influence in both theoretical and ex- 
perimental developments, proton decay, that represents its most characteristic prediction, 
has eluded all the experimental searches so far. Stringent limits on the proton lifetime in 
many channels have been set, such as [T] 

t{p 7r°e+) > 5.4 x 10^^ yr 90%C.L. 

t{p^K+v) > 2.2xl0^^yr 90%C.L. (1) 

and, as a consequence, the minimal SUSY GUTs are tightly constrained or even ruled 
out. The proton decay rate in minimal SUSY GUTs depends on many parameters and 
there are many sources of theoretical uncertainties, such as the mass of the colour triplets 
whose exchange dominates the amplitudes, the hadronic matrix elements, the spectrum of 
SUSY particles, tan f3, unknown phases and mixing angles from the flavour sector. Even 
by stretching the uncertainties to their limits, minimal SUSY SU(5) is already ruled out 
by now [2j and, in general, only non-minimal schemes, considerably more complicated 
than the minimal ones, survive this important test [H]. On the theoretical side a related 
difficulty of minimal GUTs is how to reconcile light electroweak doublets with superheavy 
colour triplets, both occurring in the same GUT multiplets: the so-called doublet-triplet 
splitting problem. In minimal GUTs such a splitting is achieved through a fine-tuning by 
fourteen orders of magnitude. Even when such a splitting is obtained, either by fine-tuning 
or by another mechanism, in many models it is upset by radiative corrections when SUSY 
is broken and/or by non-renormalizable operators which can originate from the physics at 
the cut-off scale jl]. 

An appealing mechanism to achieve the desired splitting is when the grand unified 
symmetry is broken by the compactification mechanism in models with extra spatial di- 
mensions jS]. In the last years this mechanism has been reproduced in the context of 
simple field theoretical models, GUTs defined in five or more space-time dimensions 
In the simplest case a single extra spatial dimension is compactified on a circle. The 
Lagrangian of the theory is invariant under SU(5), but the fields are periodic only up 
to transformations leaving the SU(5) algebra invariant. This requirement produces auto- 
matically a mass splitting in complete SU(5) multiplet. For instance, in the gauge vector 
boson sector, the only fields that remain massless correspond to the vector bosons of the 
standard model. All the other states, including the extra gauge bosons X of SU(5), be- 
come massive, with masses of the order of the compactification scale = 1/R. This is 
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a simple and elegant way to effectively break SU(5) down to the standard model gauge 
group. The same boundary conditions that break the grand unified symmetry should be 
consistently extended to the other fields of the theory and in particular to the multiplet 
containing the electroweak doublets and colour triplets. As a consequence, an automatic 
doublet-triplet splitting of order 1/R occurs, with no need of an ad hoc scalar sector. 

It has soon been realized that such a framework has interesting consequences also on 
the prediction of proton lifetime. Proton could even be made stable in such a construction 
[ZIIH]. In a less radical approach, the automatic doublet-triplet splitting is accompanied 
by the absence of baryon- violating dimension five operators |S] . Therefore proton decay 
can only proceed through dimension six operators, originating from the exchange of the 
superheavy gauge vector bosons X between the ordinary fermions. In conventional, four- 
dimensional models, the inverse decay rate due to dimension six operators is given by 

m 

, o+, i« Tn36 /«(GeV3)\-' (auy ( Mx (GeV) ^ 
rip Tt'e^) = 1.6 X 10 yr x — ^ x — — x ^ — — ^ , (2) 

where Mx is the mass of the gauge vector bosons X, au is the gauge coupling constant 
at the unification scale and a is a parameter coming from the evaluation of the hadronic 
matrix element. The central value in Q is too large to be observed, even by the next 
generation of experiments. In five dimensional models there are important modifications 
and the prediction, even in the most conservative case, can be quite different from the one 
in First the proton lifetime scales with the fourth power of the compactification scale 
Mc, which can be considerably smaller than the central value of Mx in eq. Q. Second, 
the coupling of the gauge vector bosons X to ordinary fermions differs, in general, from 
four-dimensional SUSY SU(5). This has to do with the freedom of introducing matter 
fields either as 'bulk' fields, depending on the extra spatial dimension, or as 'brane' fields, 
localized in a particular four-dimensional subspace of the full space-time. 

In this paper we focus on the five- dimensional SUSY SU(5) GUT described in ref. 
jni El which can be considered among the simplest realizations of GUTs in extra 
dimensions i. There are estimates of the proton decay rate in this model and its variants. 
In ref. j^l HI] the compactification scale M^ has been evaluated from the analysis of gauge 
coupling unification in a leading order approximation. By calling A the cut-off scale of 
the theory, it has been shown that, for AR equal to 100 (10), the compactification scale 
ranges in 

5 X 10^^ (3 X 10^^) GeV < < 2 x 10^^ (8 x 10^^) GeV . (3) 

Indeed AR ^ 100 is preferred by a more refined estimate Jl] and, if we naively substitute 
2Mc/'K to Mx in eq. (j21), we find a proton lifetime 

1.0 X 10^^ yr < t{j) 7r°e+) < 2.6 x 10^^ yr , (4) 

already excluded by the experimental bound of eq. ((H). In this first estimate, however, 
it has been assumed that all matter fields are coupled to the X bosons as in the four- 
dimensional theory. Actually, a more suitable setting for a natural interpretation of the 

1 Proton decay in other GUTs with extra dimensions has been analyzed in 
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fermion mass hierarchies requires that only T3 and Fj (i = 1, 2, 3) couple to the X boson 
as in the four- dimensional theory ^T] ^1 ^]. Here T and F denote 10 and 5 SU(5) 
representations and the index labels the different generations. In such a framework the 
leading baryon- violating dimension six operator is T^T^FiFi. Mixing angles are required 
to convert third generation fields into states relevant to proton decay and the amplitude 
for the dominant decay mode, K'^u, is depleted by with respect to vr^e"'" in the four- 
dimensional theory, A ~ 0.22 being a Cabibbo suppression factor. The conclusion of ref. 
[T^ is that 

t{p K+u) ^ 10^^^2 yr , (5) 

unless additional, model dependent contributions are invoked. 

The result (jSl) is not far from the expected sensitivity of future Water Cerenkov and 
Liquid Argon detectors for the K^u mode [T^ : 

0.2 X 10^^ yr lOyr, 650kton Water Cherenkov 

1.1 X 10^^ yr lOyr, lOOkton Liquid Argon (6) 

This strongly motivates a careful and complete quantitative estimate of the proton lifetime 
in this model, with a particular attention to the main sources of theoretical uncertainties, 
that can sizeably affect the estimate in eq. (0). Such an analysis is the purpose of this 
paper. First of all we compute from a next-to-leading analysis of gauge coupling 
unification, including two-loop running, threshold corrections from SUSY particles, and 
possible SU(5) breaking contributions allowed by the present construction. As we will see 
in detail in section 4, we find that there is a big theoretical uncertainty on M^, entirely 
dominated by the expected SU(5)-breaking contributions and not accounted for by the 
previous estimates. According to our analysis Mc can range from approximately 10^^ GeV 
to more than 10^^ GeV, a considerably wider interval than in eq. 0. We also size the 
uncertainty coming from the flavour sector, by contemplating different settings for matter 
fields. In what we call 'option 0', all matter fields are coupled to the X gauge bosons as in 
four-dimensional SUSY SU(5). We will also consider two other possibilities, that we call 
'option r and 'option IF. Option II corresponds to T3 and Fi {i = 1, 2, 3) coupled to X as 
in four dimensions, whereas in option I only T3 and F3 2 have standard couplings. As far as 
neutrino mass matrices are concerned, option I and option II give rise to "semianarchy" and 
"anarchy", respectively [TTj. We estimate the degree of flavour suppression by counting 
powers of the Cabibbo angle A, in two different ways, characterized by a stronger and 
a milder mixing between the different generations. This allows us to estimate also the 
theoretical uncertainty coming from the Yukawa sector. 

Our results are presented in section 5. They are qualitatively different for option and 
options I/II. In option 0, the dominant decay mode is 7T^e~^ and the lifetime is obtained 
from eq. (j21) through the substitution Mx — >■ 2Mc/7r ^. Due to the large uncertainty on 
we will see that such an option is strongly constrained, but not totally ruled out. In option 
I and II, we confirm that K~^p is one of the dominant decay channels. The prospects for 
detectability of proton decay in this channel by future machines are good, but only with 

^The numerical factor 2/tt is due to the exchange of a whole tower of Kaluza-Klein modes. 
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favorable combinations of mixing angles (option II is needed) and of SU (5)- violating terms 
at the cut-off scale. We also evaluate branching ratios for proton decay, much less affected 
by uncertainties than absolute rates. 



2 SU(5) unification in extra dimensions 
2.1 Space-time orbifold 

Following ref. [HI we consider a 5-dimensional space-time factorized into a product of the 
ordinary 4-dimensional space-time M4 and of the orbifold S^/ {Z2 x Z2), with coordinates 
x^, [fi = 0, 1,2,3) and y = x^. The orbifold 5*^/(^2 x Z2) is obtained from the circle 
of radius R with the identification provided by the two reflections: 

Z2: y^-y , (7) 

Z'^: y'^-y' iy'^y-nR/2) . (8) 

As a fundamental region of the orbifold, we can take the interval from y = to y = ttR/2. 
At the two sides of the interval, we have two four-dimensional boundaries of the space- 
time, called branes. On the covering space, S^, where we choose to work from now on, the 
origin y = and y = nR represent the same physical point and similarly for y = +ttR/2 
and y = —ttR/2. When speaking of the brane at y = 0, we actually mean the two 
four-dimensional slices at y = and y = nR, and similarly y = ttR/2 stands for both 
y = ±7iR/2. 

For a generic field (j){x'^,y) living in the 5-dimensional bulk the Z2 and Z2 parities P 
and P' are defined by 

(/)(x^,?/) ^ (f){x'',-y) = P(f){x'',y) , 

0(x^y') ^0(x^-l/') = ^'V(^^l/') • (9) 
Denoting by (f)±± the fields with (P, P') = (±, ±) we have the ?/- Fourier expansions: 



V T^R n=0 ^ 



ttR R 



^-ix'.y) = \^t<P'--^'\^n^^^^^^ ■ (10) 

ttR R 



where n is a non negative integer and the notation is such that the Fourier component field 
(p^^'^x) acquires a mass n/R upon compactification. Only has a massless component 
and only 0++ and 0^ are non-vanishing on the y = brane. The fields 0++ and |_ 
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are non- vanishing on the y — TrR/2 brane, while (j) vanishes on both branes. The 

normahzation is such that: 

/ . (11) 

2.2 Gauge and Higgs sectors 

The theory under investigation, defined in the above five-dimensional space-time, is invari- 
ant under N=l SUSY, which corresponds to N=2 in four dimensions, and under SU(5) 
gauge transformations. We can keep working in the entire circle if we also ask invari- 
ance under Z2 x Z'2 parities. The parities of the fields are assigned in such a way that 
compactification reduces N=2 to N=l SUSY and breaks SU(5) down to the SM gauge 
group SU(3) X SU(2) x U(l). Leaving aside for the moment quarks and leptons and their 
SUSY partners, the 5-dimensional theory contains the following N=2 SUSY multiplets of 
fields. First, a vector multiplet with vector bosons Am, M = 0, 1, 2, 3, 5, two bispinors A*, 
i — 1,2, and a real scalar E, each of them transforming as a 24 representation of SU(5). 



{P,P') 


field 


mass 


(+,+) 


a;, A2«, Hi, Hi 


2n 
R 


(+,-) 


A% Hi, Hi 


(2n+l) 
R 


(-,+) 


A a \1« TjT TjT 


(2n+l) 
R 


(-,-) 


Aa \la rjD tjD 


(2n+2) 
R 



Table 1: Parity assignment and masses (n > 0) of fields in the vector and Higgs super- 
multiplets. 



These fields can be arranged in an N=l vector supermultiplet V = (A^, A^) and an N=l 
chiral multiplet $ = {A^, E, A^), both transforming in the adjoint representation of SU(5). 
Then there are two hypermultiplets H^ (s — 1,2), each of them consisting of two N=l 
chiral multiplets: H^ — {H^,Hi) and H'^ — {Hi,H^). Here Hr, and Hi stand for the 
usual SU(5) supermultiplets, which include the scalar Higgs doublets Hi and Hi and 
the corresponding scalar triplets Hi and i?J, respectively. The other two, H^ and H^, 
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transform as 5 and 5, under SU(5). The parity P, P' assignments are the same as in ref. 
jn] and are given in table 1. 

Here the index a (a) labels the unbroken (broken) SU(5) generators T°- (T°), H stands 
for the whole chiral multiplet of given quantum numbers. The parity P breaks N=2 SUSY 
down to N=l and would allow complete SU(5) massless supermultiplets, contained in the 
first two rows of table 1. The additional parity P' respects the surviving N=l SUSY 
and breaks SU(5) down to the standard model gauge group. Note that the derivative 

transforms as (— , — ). The (+, +) fields, which remain massless and do not vanish on 
both branes are the gauge and Higgs multiplets of the low energy Minimal SUSY Model 
(MSSM). The bulk 5-dimensional Lagrangian is exactly as in ref. jOj and we do not 
reproduce it here. 

The vector bosons ^4^^°^ remain massless and, together with the gauginos A^"^'^-', form 
a vector supermultiplet of N=l. All the other vector bosons are massive and become 
component of N=l massive vector supermultiplets. Therefore the SU(5) gauge symmetry 
is effectively broken down to SU(3) x SU(2) x U(l) and the symmetry breaking order 
parameter is the inverse radius 1/-R, which is expected to be of the order of the unification 
scale. At the same time, in the Higgs sector, only the Higgs doublets and their super- 
partners are massless, while colour triplets and extra states acquire masses of order 1/-R, 
giving rise to an automatic doublet-triplet splitting. The N=2 SUSY algebra would allow 
a mass term between the hypermultiplets and thus spoiling the lightness of 

the Higgs doublets achieved by compactification, but such a term can be forbidden by 
explicitly requiring an additional U(l)ii; symmetry [HE]. This U(l)/j, symmetry is the 
diagonal subgroup of an i?-symmetry acting on the Higgs multiplets as 

{H,,H-,)'{e) = e^-(i75,i^5)(e-^"^) (12) 
{H-,,H,)'{e) = e^^{H-,,H,){e-^'^e) (13) 

and a U(l) symmetry acting as: 

{H,,Hl)' = e-^%H,,Hl) (14) 
{H-,,Hl)' = e-"^{H-„Hl) (15) 

where 6 denotes the Grassmann coordinate of the N=l superspace. The overall R charges 
of Hq, ifg, ifg and are 0, 2, 0, 2 respectively. This symmetry forbids a mass term 
between Hr, and H^. Therefore, before the breaking of the residual N=l SUSY, the mass 
spectrum is the one shown in table 1. 



2.3 Matter fields 

Different options exist to introduce quarks and leptons. A first possibility is to describe 
them as N=2 bulk hypermultiplets jB]- In each hypermultiplet the different SU(5) compo- 
nents undergo a sphtting similar to the one occurring in the gauge and in the Higgs sector. 
Consider for instance a hypermultiplet transforming as 5 © 5 under SU(5), where, in an 
obvious notation, 5 = {L,D'^) and 5 = {L,D^). We can assign orbifold parities in such a 
way that either L or D'^ but not both of them contains a fermion zero mode and, to get all 
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the massless degrees of freedom of a 5 representation, we need to introduce an additional 
hj^ermultiplet, (L, D'^)' Q) (L, D'^)' with opposite P' assignment. A similar splitting occurs 
in 10 representations. It is also interesting to note that in the bulk Lagrangian, that con- 
serves the momentum along the fifth dimension, the gauge vector bosons that violate 
baryon number, never couple two fermion zero modes of an hypermultiplet. Therefore 
fermions belonging to hypermultiplets cannot contribute to proton decay, at least in a 
minimal version of the model ^. 

A second possibility is to identify matter as N=l SUSY chiral multiplets, localized at 
one of the two branes. On the brane a.ty = ttR/2, the gauge invariance is limited to the SM 
gauge transformations and incomplete SU(5) multiplets can be introduced. The kinetic 
terms of such multiplets, localized at y = ttR/2, do not contain any baryon or lepton 
violating interactions. Therefore such interactions are only possible if we introduce ad 
hoc non-minimal terms of the gauge vector bosons in y = ttR/2. Although technically 
possible, we will disregard such a possibility in the present analysis, not to introduce too 
many unknown parameters. On the brane at ?/ = 0, where the full SU(5) gauge invariance 
is at work, N=l chiral multiplets should be introduced as complete SU(5) representations. 
Moreover, the gauge invariant kinetic terms defined on the brane do not conserve the 
momentum along the extra dimensions, and the fermions are coupled to all the gauge 
vector bosons, including those violating baryon and lepton numbers. In summary, in 
a minimal model, matter is introduced either in bulk hypermultiplets or in N=l chiral 
multiplets at ?/ = 0. Baryon and lepton violation between massless fermions mediated by 
gauge vector bosons is only due to minimal couplings at y = 0. 

Additional sources of baryon and lepton violations can be provided by Yukawa cou- 
plings that, in the present model, are all localized in y = 0. Among these couplings there 
are dimension three and four operators, which in MSSM are forbidden by R-parity. The 
role of R-parity is played in the present context by the R-symmetry discussed in 

the previous subsection, forbidding a mass term between the Higgs hypermultiplets 
and H^. As shown in ref. [HI and described in table 2, such a symmetry can be suit- 
ably extended to the matter sector in such a way that the leading order allowed Yukawa 
interactions are: 

Sy = j d'xdy5{y) j d^9 (ry^TH, + Fy.TH-, + ^^']^^^' ) + h.c. (16) 

Here T and F denote 10 and 5 representations, either localized in y = or belonging to 
bulk hypermultiplets and yu^d and w are the Yukawa couplings, 3x3 matrices in flavour 
space; M denotes the scale where the total lepton number is violated. In the limit of exact 
residual N=l SUSY, the symmetry is also exact and all the unwanted dimension 

3 and 4 operators are absent. When N=l SUSY is eventually broken, such operators are 
still forbidden if R-parity, contained in U(l)ij, survives after the breaking. 

An important consequence for matter stability is that the dimension 5 operator TTTF 
is also forbidden by U(l)/j. Indeed, a crucial difference between this model and minimal 

■^It is possible to introduce non-diagonal kinetic terms between matter hypermutiplets (M, M') on 
the y = brane ^]. Such terms can contain interaction terms between matter zero modes and 
superheavy gauge bosons, thus contributing to proton decay. In our analysis we neglect such a possibility 
that enhances proton decay rates at the price of introducing a strong model dependence. 
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Field 


V E 


H5 


H5 


T T 


F F 


U{1)r 





2 


2 


1 1 


1 1 



Table 2: U(l)/{ charges for 4D vector and cliiral superfields. Matter fields, either in 
tenplets T or fiveplets F of SU(5) have equal U(l)i^ charge, whether they are bulk or 
brane multiplets. 

SUSY SU(5) is that the colour triplets whose exchange gives rise to TTTF acquire mass 
from compactification and not from a term that mixes and ifg. When U(l)/j gets 
broken at low energy a mixing term of the order of the SUSY breaking scale arises, but 
it is too tiny to produce significant effects. As a result, the only sizeable source of proton 
decay in the model are the dimension 6 operators originating from the exchange of the 
heavy fields A^. 

The Yukawa interactions of eq. (fTBjl describe, after electroweak symmetry breaking, 
masses for fermions, including neutrinos. A careful inspection of fermion masses suggests 
whether to localize matter multiplets at ?/ = or to embed them in bulk hypermultiplets. 
The most plausible assignments are reported in table 3. 





option 


option 1 


option 11 


bulk 








brane y = 


T„ Fi (z = 1,2,3) 


^3, F2, F3 


T3, Fi, F2, F3 



Table 3: Matter fields and their locations. Indices 1, 2, 3 label generations. Bulk fields 
belong to N=2 supermultiplets and they should be doubled, to provide the correct number 
of zero modes. In the table this doubling is understood. For instance Ti of option I stands 
for Ti, T{. Options I and II differ only for the assignments of Fi. 



We shortly recall the origin of this choice. Option realizes the simplest five-dimensional 
realization of SUSY SU(5). The hierarchy among fermion masses is accommodated by ad- 
justing by hand the Yukawa couplings. The same fermion mass relations of minimal SUSY 
SU(5) hold. Option I and II represent alternative description of the flavour sector that 
take advantage of the existing freedom in five dimensions. If both Ti and Fi of the same 
generation i are put in y = 0, then, as in minimal SUSY SU(5), this implies the equality 
between the ii mass matrix elements of charged fermions and down quarks. This is only 
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acceptable, and actually welcome, for the third generation. Hence, to begin with, in option 
I and II, T3 and are chosen as brane fields. On the contrary, if at least one between T 
and F belong to an hypermultiplet, the SU(5) relation does not hold any more, because 
the masses of charged fermions and down quarks originate from two different hypermulti- 
plets. By assigning F2 to y = we can naturally explain the largeness of the 23 mixing in 
the lepton sector, and, consequently, T2 is chosen as a bulk field. It is also preferable to 
have Ti as a bulk field, which helps in reproducing the strong hierarchy in the up quark 
sector. The only ambiguity concerns Fi and we are left with the two possible options, I 
and II, of table 3. 



2.4 Fermion masses 

Fermion masses in option are introduced by tuning the Yukawa couplings ?/„, yd and w 
of eq. (fTBj) . to the desired values. The resulting description is not completely satisfactory 
due to the equality between down quark and charged lepton masses, the well known mass 
relations of minimal SU(5). These are good for the third generation but wrong, by factors 
of order one, for first and second generations. There are several mechanisms to correct 
the undesired relations [19^ but they do not appreciably alter the predictions for proton 
lifetime [20] • Therefore, in what follows, we will maintain option in its minimal form, 
without introducing any Yukawa couplings other than those of eq. flTBj) . 

Alternatively, in option I and II, we can take advantage of the fact that bulk and brane 
matter fields enter the Yukawa couplings of eq. with a different relative normalization 
which is induced by the different mass dimension of four and five dimensional fields. Each 
Fourier component of a bulk field in eq. ^TE\i carries a factor 

^ (17) 



as can be seen from the expansion of eq. ()1()|1 . The mass parameter A, needed for dimen- 
sional reason in the Yukawa couplings, is some cut-off scale of our effective theory. Notice 
that e is always a suppression factor since the theory only makes sense as an effective the- 
ory in the regime A > 1/R. We can include the appropriate power of e into the coupling 
constants yu,d and w. From the assignments of table 3 we can read the relative suppression 
due to e of the different matrix elements in fermion mass matrices: 



Vu 









e 




1 J 
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option I option II 





( 


e2 


e ^ 




( e 


e 


1 


\ 


C/LL 


e 


e 


1 




t 


e 


1 








e 


1 , 




V e 


e 


1 


/ 

/ 






e 






^ 1 


1 


1 


\ 


W = 


e 


1 


1 


w = 


1 


1 


1 








1 








1 


1 





where the matrices are given up to overall factors and up to further order-one dimension- 
less coefficients that might affect each individual entry. The charged lepton mass matrix 
m-e is approximately given by mj, the exact equality being broken by the first and second 
generations. An interesting structure shows up in the 23 sector. If e is of order A^, where 
A ~ 0.22 denotes a Cabibbo suppression factor, then mc/nit ~ A^, nis/mh A^, and 
Vcb ~ A^. In the lepton sector we have m^/rrir ~ A^, ^23 ~ 0{1) and all this goes into 
the right direction. From eq. ()18I19|) we see that the first generation is not equally well 
described. We can adjust first generation couplings by exploiting the freedom related to 
the dimensionless coefficients, which in eq. (118119^ have been set to one. Alternatively we 
can improve the description of fermion masses by slightly modifying the wave function for 
the zero mode of Ti with respect to that of T2. So far, these wave functions have been 
taken constant in y with the same value at ?/ = 0. If the wave function for the zero mode 
of Ti is displaced from y = there will be a relative suppression between the couplings of 
Ti and T2. Such a change in the ^/-profile of the wave function for the zero mode in Ti can 
be easily obtained by introducing a bulk mass term for the hypermultiplet containing Ti 
and the net effect can be described by replacing e with a new suppression factor 5 < e in 
the matrix elements involving Ti [21]. We get: 



option I 



Vu 



/ 52 e5 5\ 
6 e 1 



(20) 



option II 



f eS e2 e \ 



Vd 



6 
6 



e 1 
e 1 



Vd 



fS e 1\ 
6 t 1 
5 e 1 



(21) 



By assuming that the unknown coefficients multiplying the different matrix elements are 
of order one, we get the order of magnitude relations: 



Tfl^. 2 r-2 

— K, e — ^ 

rrit rrit 

Ills rrid rrie 
nil, [ S (optionll) 



e6 (optioni) 

(22) 
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s 



ub 



cb 



(23) 



We can match almost all experimental data, by taking e of order and 6 of order A'^ -i- A'^. 
We observe that the choice S = seems more appropriate to reproduce the observed 
quark mixing angles, needed for our estimate of the proton lifetime. For neutrinos, in 
both options I and II, we should invoke some additional accidental relations not 
automatically guaranteed by the textures in eq. ()19|). In case I, known as "semianarchy" , 
the determinant of the 23 block in m^, which is generically of order one, should be tuned 
around 1122/1713 ^ J Aml^J Aml^^ ^ O.l-i-0.2. We should also enhance somewhat the solar 



mixing angle. In case II, denoted as "anarchy" [25: , we need again the same cancellation in 
the 23 determinant and we should adequately suppress the 6^13 angle. The overall picture is 
not entirely satisfactory, though, at least from a technical view point, the modest required 
tuning can be easily realized by exploiting the freedom in the dimensionless coefficients 
that multiply each matrix element in the above textures. 

When computing the amplitudes for proton decay, we will need the transformations 
mapping fermions from the interaction basis into the mass eigenstate basis. We can express 
these transformations as: 



u 



u 



LuU 



d 



Ldd 



Lee 



(24) 



where in the left-hand (right-hand) side fields are in the interaction (mass eigenstate) basis 
and L, R are 3x3 unitary matrices. The quark mixing matrix Vckm and the lepton one 
UpMNS are given by: 



Vr 



CRM 



U 



PMNS 



LlL. 



The unitary matrices L and R can be estimated from eqs. 
(Uni). We find: 



(25) 

and the action in eq. 



Lu — Ru 



Rp. 



f 1 X f 1 X^ X^\ 

X 1 A2 
\ A=^ A^ 1 J 



X^ 1 A2 
V A^ A^ 1 J 



(26) 



letting S vary in the range A^ -i- A^. When 5 = A^ some accidental enhancement is required 
to correctly reproduce Vus and Vub- We have also 



option I 



Le ~ Rd 



( 

X^ 



1 A2 A2 \ 

1 1 



option II 



1 1 



Rd ~ L,, 



T ~ 



I I X X^\ 
X 1 1 



/l 1 

1 1 1 
VI 1 1 



1 \ 



(27) 



1 1 
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where, in option I, we have assumed that the determinant of the 23 block of w is of order 
A. Moreover, in option II, some accidental cancellation should guarantee 6*13 ~ A. Eqs. 
()26I27|) are basic inputs in our estimates of the proton lifetime. 

2.5 Supersymmetry breaking 

In the model under discussion, baryon violating processes are dominated by dimension 
six operators originating from vector boson exchange and, in first approximation, the 
proton lifetime is not sensitive to the breaking of the residual four- dimensional N=l SUSY. 
Nevertheless a dependence on the light supersymmetric spectrum arises from the threshold 
corrections to gauge coupling unification, around the SUSY breaking scale. As we shall see 
in detail in section 4, the compactification scale that controls proton decay amplitudes is 
determined by gauge coupling unification and, in a phenomenological analysis, is slightly 
affected by the light supersymmetric particles. 
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Table 4: Spectrum of SUSY particles. Masses are in GeV. The values are taken from ref. 
and refer to different realistic scenarios. Here they are given in an SU(2) invariant approximation. 



There are several SUSY breaking mechanisms that can be adapted to the present setup. 
It is possible to break SUSY by non-trivial boundary conditions on the bulk superfields 
122 We will denote such possibility as Scherk-Schwarz SUSY breaking (SSSB). In 
this case, at the tree-level, the splitting within N=l supermultiplets only affects bulk 
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610 


425 


563 


214 
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Table 5: Masses in GeV of SUSY particles within SSSB, from ref. in a SU{2) invariant 

approximation. In the sfermion sector the first row refers to first and second generation, while 
the second row corresponds to the third generation. 



fields and is proportional to 1/R, which in the present context is close to the unification 
scale. Therefore, to keep superparticles at the TeV scale, it is necessary to have a tiny 
dimensionless SUSY breaking parameter [23 ■ SUSY can also be broken by an intrinsic 
four-dimensional mechanism on either of the two branes. The SUSY breaking sector and 
the transmission of SUSY breaking to the observable sector are notoriously a source of 
ambiguities and phenomenological problems. Here we are only interested in the estimate of 
the theoretical uncertainty on proton lifetime coming from our ignorance about the SUSY 
breaking mechanism. We think that, from this viewpoint, the best way of parametrizing 
such an ignorance is to assume a variety of supersymmetric particle spectra and, for each 
of them, evaluate the proton decay rate for the various channels. The comparison among 
the results corresponding to different spectra will quantify the theoretical uncertainty we 
are interested in. As we will see, this uncertainty is of a certain relevance but it is not the 
dominant one in our problem. 

In our numerical estimates we have adopted the ten spectra listed in table 4. These 
spectra have been adapted from ref. [211 where they correspond to the so-called Snowmass 
Points and Slopes (SPS), a set of benchmark points and parameter lines in the MSSM 
parameter space corresponding to different scenarios. We have simplified the spectra of 
ref. j^ni, by taking an SU(2) limit where all left-handed sfermions are degenerate. Similarly 
winos and higgsinos are taken as unmixed mass eigenstates. In addition, in our option II 
for matter fields, we have also considered the spectrum of ref. ^2); obtained by breaking 
SUSY through boundary conditions for the bulk fields. This spectrum is listed in table 5. 



3 Dimension 6 operators for proton decay 

The interactions contributing to proton decay are those between the gauge vector bosons 
X = A'J^ and the N=l chiral multiplets on the brane at y = 0. In option 0, where all 
matter fields are at y = 0, the effective lagrangian for proton decay is formally identical 
to that from gauge boson exchange in minimal SUSY SU(5): 



^ 2M2. 
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- ul'j''Ui ■ dp^Cj + u1'-f^{VcKM)ijdj ■ dl'j^Uk , (28) 

where all fermions are mass eigenstates (neutrinos are taken as massless) and i,j,k = 1,2,3 
are generation indices. Colour indices are understood. The dimensionless quantity 

is the four-dimensional gauge coupling of the gauge vector bosons zero modes. The com- 
bination 

M.^?^ . (30) 

TT 

proportional to the compactification scale 

Mc^^ , (31) 

is an effective gauge vector boson mass arising from the sum over all the Kaluza-Klein 
levels: 

°° 1 1 

S {2n + lyM^ = 2M| • 

When expressing interaction eigenstates in terms of mass eigenstates, through eq. (j21I), 
we have exploited the minimal SU(5) relations Kl = L„, i?^ = Lg, R\ = La, which hold 
also in option I. 

When considering option I and II the interactions between and the matter multiplets 
on the brane aX y = are 

Sab^o = / c?'a;4+-)"^(x, 0) J^,(x) + h.c. , (33) 

where A^^ ^"\x,y) is the five dimensional gauge vector boson with quantum numbers 
(3,2,5/6) under SU(3)x SU(2) x U(l), a and / are color and SU(2) indices respectively. 
The currents Jj^j are given by: 

Jal = (-aM <^Ijq1^ Y + e^3 Y <lan " d'^aa 7^ ha 

= q3l''ul+^3Yq3-dr-aYla , (34) 

where in the second line gauge indices have been abolished and only generation labels are 
presents: the sum over a is restricted to the multiplets Fa living on the brane at y = 0. 
All fermions q, I, u^, d^ and e'^ are left-handed in our notation. If we insert the expansion 
of eq. (fTUI) and we integrate out the super heavy gauge bosons in the limit of vanishing 
momenta for the light particles, we obtain the four-fermion Lagrangian: 

,2 



= {uirqz ■ el^.q, - u^^qz ■ d%i,la) • (35) 

By expressing it in terms of mass eigenstates through the transformations of eq. ()24p , we 
get: 

2 



= -jgr {^1 «'7^w ■ et^^d, - Cl u-^^u ■ d]^^e, + Cj, u-^^di ■ dp^u^) . (36) 
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Here u denotes the lightest up type quark, i,j = 1,2 and k = 1,2, 3. The coefficients C"s 
are given by: 



where, as above, the sum over a is restricted to the multiplets Fa hving on the brane at 
y = 0. From the previous equations we can easily estimate the leading contribution to the 
decay amplitudes. In option the result coincides with those due to X boson exchange in 
minimal SUSY SU(5), with the substitution Mx — ^ 2Mc/tt. In particular the dominant 
decay mode is 7T^e~^, whose branching ratio is approximately 30 -r- 40% fI7\ . 

In option I and II, by neglecting any numerical coefficients and any subleading contri- 
butions, we obtain the results listed in table 6. From table 6 we see that the contribution 
from the first term in the expression of the Lagrangian proportional to the coeffi- 

cients Clj, is always negligible. This is due to the high suppression brought by the mixing 
matrices when converting third generation fields into light fermions. From this point of 
view the suppression is less severe for the second term in eq. ()35|) . which gives origin 
to Cfj and Cf^f^. When Fi is taken as a bulk field (option I), the dominant operators are 
uusfi and udsv^^r-, with a fiavour suppression of order A^^^. This favors the decay channels 
p —>■ K^fi~^ and p —>■ K^v. When Fi is chosen as a brane field (option II), the dominant 
operator is udsv, with a neutrino of whatever type and a suppression of order A^^^. The 
preferred decay channel is K^v. All the other possible operators contributing to proton 
decay, listed in table 6, are depleted with respect to the dominant one, with a parametric 
suppression A^^^. This discussion shows that, even in the most favorable case, we should 
expect a large fiavor suppression in the decay amplitude, with respect to option where 
essentially no small angles are required to proceed from fiavour to the relevant mass eigen- 
states. Numerically A^^^ ~ (1 ^ 5) x 10~^ and it would seem quite difficult to rescue such 
a small numerical factor and end up with a measurable effect. In the next section we will 
carefully analyze gauge coupling unification in this class of models, in order to quantify 
the overall coefficient gij/{2Mx) and its impact on the proton lifetime. 

4 Gauge coupling unification 

The overall strength of proton decay amplitudes is governed by the masses of the gauge 
vector bosons A^, which, in the present model is proportional to the compactification scale 
Mc- The latter can be estimated by analyzing the unification of gauge couplings. The 
low-energy coupling constants aj(m^) {i = 1,2,3) in the MS scheme are related to the 
unification scale Au, the common value au = gfj/^Aix) at A^/ and the compactification 
scale Mc by the renormalization group equations: 




2{Ru)u{Lu)3l{Re)i?,{Ld)3j 




(37) 
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option II 
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Cf, ^ A6-8 








uuse 






UUSyU 




Cl ^ A6-« 


uddz/g 

u.u.u.iy^ 7- 


Cf,, ^ A^-ii 

{k=2,3) 


I ft, — J_ j^^O J 


udsi/g 


(A;=2,3) 


CI,, ^ A^-6 

(fc=l,2,3) 



Table 6: Leading contributions to the different decay amplitudes from the Lagrangian of 
eq. (jnni)- Only the order of magnitude in A is displayed. Powers of A vary in the range 
displayed in the table depending on the choice 6 ^ -r- X'^. 
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Here bi are the coefficient of the SUSY /? functions in the one-loop approximation: 



(b, ] 




( 33/5 \ 


b2 




1 


V ^3 ) 




I -3 ) 



(39) 



for 3 generations and 2 hght Higgs SU(2) doublets. We recall that gi is related to the 
hypercharge coupling constant gy by gi = ^5/3 gy- In eq. flHHj) . 6i stand for non-leading 
contributions and depend upon Mc- More precisely: 



6, 



('2') L 

where 61 represent two-loop running effects, 61'' are light threshold corrections at the 
SUSY breaking scale, 6^^^ are heavy threshold corrections at the compactification scale 
Mc and 6f''' denote SU(5) breaking contributions from the brane ai y = tcR/2. By setting 
5j = we go back to the leading order approximation where the three gauge coupling 
constants exactly unify at the scale Au- In such an approximation, from the experimental 
values of aemi^nz) and siia^ 9w{mz) we can derive Au, au and as^mz). In the present 
context, we will use equations (jnH|) to evaluate contained in 6l'^\ starting from the 
experimental values aemi'fnz), sin^ Owi'mz) and 0:3(771^) and from estimates of ^f^'*, 6f^ 
and We use the experimental inputs 



S?^ + + ^ + 6? 
(0 



(40) 



sin^ Owi^z) 
as^mz) 



127.906 ±0.019 
0.2312 ±0.0002 
0.1187 ±0.0020 



(41) 



4.1 Two-loop 

We include the two-loop corrections coming from the gauge sector 



s(2) 



^ i=i b, 



b, 



'3 — 8 sin^ 9w 
36 sin^ 6w — 3 



(42) 



where the matrix b^j is given by: 



( 199/100 27/20 22/5 \ 

9/20 25/4 6 
^ 11/20 9/4 7/2 y 



(43) 



4.2 Light thresholds 

The threshold effects from the SUSY breaking scale are provided by j30j : 



A« = -ij:6f(,)iog(^) 



(44) 



where the index j runs over the spectrum of SUSY particles of masses mj and extra Higgses 
and the coefficients bf\j) are given in table 7. We work in the approximation where SU(2) 



17 



breaking effects are neglected. The meaning of bf\j) is very simple and can most easily 
be captured in the approximation where all particles have a common mass msusY- In 
such a case, from table 7 and eq. (j44j) we obtain: 

r{i) 5 rnsusY m) 25 msusY Ai) 2 msusv (.r-\ 

5\> = -—\og 5^ = - log 5l> = — log (45) 

47r mz 127r mz vr mz 
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1/12 




1/30 





1/12 
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1/20 


1/12 







1/20 


1/12 





extra H 


1/20 


1/12 






Table 7. Coefficients hi (j) for SUSY particles and extra Higgses. All SU(2) breaking 
effects are neglected. Sfermion contributions in the table refer to one generation. 

If we add these contributions to the leading order terms, we get 

1 1 ^ hf^ , ( msusY \ ^ h , f Au \ i.^. 

—( 7 = — + ^ log + , 46 

OLi[mz) au 27r \ mz J 27r KmsusyJ 
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where {bf^', 6f^, b§^) = (41/10, -19/6, -7) are the coefficients of the f3 functions in the 
one-loop approximation in the standard model. In this simplified case the running of the 
gauge coupling constants is standard from mz to msusY and supersymmetric from uisusy 
up to the unification scale Ku- As explained in section 2.5, we do not adopt a specific 
scheme for SUSY breaking. In the model under discussion, where dimension 5 baryon- 
violating operators are forbidden, the coefficients 5f^ are the only quantities entering the 
evaluation of the proton lifetime that are sensitive to the details of SUSY breaking. To 
the purpose of providing a realistic estimate of the proton lifetime, we can better size 
the uncertainty from the SUSY breaking sector by considering the variety of sparticle 
spectra listed in tables 4 and 5, coming from several mechanisms. As a consequence, the 
contribution (5f ■* varies in a range that describes our ignorance about the SUSY breaking 
sector. 



4.3 Brane kinetic terms 

The contributions df* in eq. (j4(jp come from kinetic energy terms for the gauge bosons of 
SU(3), SU(2) and U(l) on the brane at ?/ = 7r_R/2. These terms, which break SU(5), are 
allowed by the symmetries of the theory and, even if we set them to zero at the tree-level, 
they are generated by radiative corrections 1^ . The net effect of these terms is to 
modify the boundary value of the gauge coupling constants (?i(A) at the cut-off scale A: 

1 _ 2ttR 1 

where (75 (A) is the SU(5)-invariant gauge coupling constant, having mass dimension - 
1/2, coming from the bulk kinetic terms of the gauge bosons, while g(^i,)i are dimensionless 
gauge couplings arising from independent kinetic terms of the SM gauge fields on the SU(5) 
breaking brane at y = ttR/2. If the SU(5) breaking terms l/g'^^^ where similar in size to 
the symmetric one, we would loose any predictability and we could not hope to reasonably 
constrain the compactification scale Mc, through the evolution of gauge couplings. It has 
been observed that a predictive framework is recovered if we assume that the theory is 
strongly coupled at the cut-off scale A[/ [TT]. In such a case, dimensional analysis suggests 
that gl{A) ^ 167r^/A, gf^^^^ ^ lOvr^ and eq. (jUj) becomes: 

' , (48) 



^2 (A) 87r2 'I67r2 

Since Ai? ^ 1 is expected for the consistency of the model and Ai? ^ O(IOO) is welcome 
from the view point of fermion masses, we see that the SU(5) invariant component indeed 
dominates over the non-symmetric one and that the theory predicts gauge couplings of 
order one at the cut-off scale. In what follows we will adopt the assumption of strong 
gauge coupling at A. Of course, since our quantitative analysis will allow to exphcitly 
evaluate both R and A, we will be able to check the validity of this assumption. For our 
numerical analysis, we will regard df* as random numbers with a flat distribution in the 
interval (-l/27r, -hl/27r): 



5f^ e 
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Given our ignorance about the ultraviolet completion of our model, where we could predict 
the parameters df'^ , each value in the above interval has the same "probability" . This will 
induce an uncertainty in the determination of the compactification scale Mc, which, as we 
shall see, represents the largest theoretical error. 



4.4 Heavy Thresholds 

Two different approaches to the evaluation of the heavy threshold effects exist in the 
literature. In ref. [H], which we will refer to as HN, these effects are evaluated in a leading 
logarithmic approximation for the particles whose masses are smaller than the cut-off scale 
A. These states belong to the Kaluza-Klein towers of the gauge bosons and the matter 
hypermultiplets. The result is: 

N \ N \ 

cih)HN 7i 1 , 1 r^n\ 

= 2? S wTm. ^ S (2;7tt)m: 

where the coefficients 7^ and r^j are listed in table 8. 

The sums stop at such that (2A^ + 2)Mc is the Kaluza-Klein level closest to, but 
still smaller than, the cut-off A: 

(2iV + 2)^A . (51) 
It is convenient to recast eq. (jSUj) in the form 



^r" = 4'^^ + ^Elog#^ , (52) 



where Oi are given in table 8 and 5[f ^ is a universal contribution that can be absorbed into 

OLV'- 

^{j^) = _ (-4 + ng + 3nio) y log — r , (53) 

^ 27r^ ^ ' ^ (2n-f- l)(2n + 2)M2 ' ^ ' 

where ng and nio are the number of matter hypermultiplets transforming as 5 and 10 of 
SU(5). For large A^, that is for Ai? > 1 

Y bg -, ^ — bg A^ + 1 - - bg TT ^ — bg bg - . 54 

^0 ^ 2r2 + 2 2 ^ 2 ^ 2 ^ Me 2 ^ 2 ^ ' 



-j^log- + ... , (55) 



In this limit: 



where dots stand for universal contributions. From eq. (j50l52j) we also see that, up to an 
irrelevant universal contribution redefining the initial condition «[/, all the effect comes 
from the shift between even and odd Kaluza-Klein levels that removes the degeneracy 
within full SU(5) multiplets. 

In ref. ^3], referred to as CPRT, df^^ are estimated within an effective lagrangian ap- 
proach. The full towers of Kaluza-Klein states are integrated out in the higher-dimensional 
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model and an effective theory at the scale is obtained. The gauge couplings at the M^. 
are related to the gauge couplings at the cut-off scale by the threshold corrections: 

5^ = log - 1 - ^) - ^ log2 + - (56) 



where 



1 r+oo , „ 

1=-l dt{t-^ + e^"*" - 1) ^ 0.02 , (57) 



and 5^ are conversion factors from the DR to the MS scheme |32j also given in table 8. 

In the limit A 3> Mc, the term proportional to log(A/Mc) dominates both ^9^^^^^ in 
eq. and ^^^^^^^^'^ in eq. (j^. Actually, since 7^ = — 4, up to an irrelevant universal 
contribution, the two expressions of eqs. and ()56p coincide in this limit. In our anal- 
ysis we have used both ^f'^^'^^ and ^^'^^^^^'^ to estimate the compactification scale Mc. 
The difference between the values of obtained in these two approaches is part of the 
theoretical error affecting our estimates. 
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Table 8. Coefficients for the heavy threshold corrections due to the gauge and Higgs 
sectors eqs. (I50l52l56j) . from ref. and [Hj. If rig and Uiq matter hypermultiplets 
transforming as 5 and 10 of SU(5) are present, the additional contribution (ng -|- 3nio) 
should be added to both 7, and rji, for all i = 1,2, 3. (Two copies of n SU(5) multiplets 
give rise to n zero modes in the present model.) 

We can have a qualitative understanding of the non-leading effects by looking at the strong 
coupling constant 03 (m^). At leading order, defined by setting 6i to zero in eq. from 
the input values of aem('^z) and sin^ 9w{^z) we obtain the prediction: 



a 



LO 



(mz) ^ 0.118 



(5^ 



in excellent agreement with the experimental value. When we turn on the non-leading 
corrections 6i, the predicted value of a3{mz) is changed into: 



5. 



a 



7 



LO 



{mz) 



a 



LO 



{mz)Ss 



(55i - 12^2 + 753 



(59) 
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where the combination 6s is the sum of the two-loop correction, 6^ , the hght threshold 
correction, 5^^^ , the brane contribution and S^'^^ , the heavy threshold correction. From our 
previous discussion, these contributions are readily evaluated and we obtain: 



si'' 



-0.82 



0.50 + 



287r 
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log 



T^SUSY 



niz 



1 




(60) 



where we have considered a 'typical' spectrum of supersymmetric particles and only the 
leading logarithmic correction from heavy thresholds. The combined effect of two-loop 
corrections and light thresholds would raise the prediction of 03 (m^) up to approximately 
0.129, exceedingly larger than the measured value. This can be brought back inside the 
experimental interval by the corrections from heavy thresholds, provided there is a gap 
A/Mc ~ 100. This qualitative discussion shows that in the present model gauge coupling 
unification can take place with a compactification scale considerably smaller than the cut- 
off of the theory and, hopefully, smaller than the unification scale in conventional GUTs. 



We have performed a numerical analysis in order to evaluate the compactification scale 
Mc- Inputs to our computations are the experimental results in eq. (PT|) and the spectra of 
supersymmetric particles and extra Higgses given in section 2, tables 4 and 5. The brane 
contribution sf'^ described in section 4.3 are taken as random variables in the interval of eq. 
()49|) . where they are generated with a flat distribution. Fig. 1 shows the compactification 
scale Mc as a function of the spectrum of supersymmetric particles. The heavy threshold 
corrections are evaluated according to both the approaches HN (error bars on the right) and 
CPRT (error bars on the left) discussed in section 4.4. For each given SUSY spectrum there 
are two main sources of errors. One is the experimental error, approximately Gaussian, due 
to the uncertainties affecting the input parameters in eq. ()41|) . This is largely dominated by 
the error on 0:3(771^) and is represented by the smallest error bar in the figures. The other 
one is theoretical, associated to the unknown SU(5)-violating kinetic terms aX y = 7ri?/2. 
This is a non-Gaussian error, since we have no reason to prefer any values of the parameters 
in the interval (jl^ . The linear sum of these two errors is described by the largest 
error bar in our figures. The total error is therefore fully dominated by the theoretical 
one and it heavily affects the prediction of the compactification scale, that is predicted 
only up to an overall factor 10^^. The prediction of Mc also depends on the spectrum of 
SUSY particles at the TeV scale. The induced error on Mc is however subdominant and 
the largest difference between the values of Mc obtained by varying the SUSY spectrum 
is by a factor of 4, at most. Finally, also the treatment of the heavy thresholds gives 
rise to a theoretical uncertainty. The values of Mc obtained by the HN procedure are 
systematically smaller, by approximately a factor 2-=- 3, with respect to those given by the 
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Figure 1: Compactification scale Mc versus SUSY spectrum. For each SUSY spectrum, the 
bigger (smaller) circle on the right (left) shows the result in the HN {CPRT) approach. The 
shorter error bar represent the parametric error dominated by the experimental uncertainty 
on aslmz), the wider bar includes the dominant source of error, the SU(5)-breaking brane 
terms. 
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Figure 2: Gauge coupling au = gfj/^An) at the compactification scale versus SUSY spec- 
trum. Error bars as in fig. ^ 

CPRT approach. Despite the large overall uncertainty, from fig. 1 we see that on average 
the compactification scale is smaller than the unification scale of four-dimensional 
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SUSY SU(5), ^ 2 X 10^^ GeV, by a factor of order 10. The effective gauge boson mass, 
Mx = 2Mi./tt, entering tfie dimension six operator of eqs. ()28I36|) is smaller than the 
corresponding mass in four- dimensional SUSY SU(5) by a factor 15. This produces 
an average enhancement by a factor 5x10^ in the proton decay rate. Moreover, the key 
point in the present analysis is that the theoretical uncertainties are completely dominated 
by the unknown SU(5) violating brane interactions and by the SUSY spectrum. Given 
the present knowledge, we cannot prefer any brane interactions or any SUSY spectrum 
among the various possibilities and the average prediction has not the meaning of the most 
probable one. From this viewpoint values of as small as 10^^ GeV are equally probable 
than Mc ~ 10^^ GeV, and the enhancement of the proton decay rate can be as large as 
5x10^. Such an enhancement is sufficient to overcome the huge suppression factor coming 

5 I 1 
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Figure 3: Ratio KjMf. versus SUSY spectrum. Error bars as in fig. ^ 

from fiavor mixing in options I and II, which in the most favorable case, as we can see 
from table 6, is |A^p 3 x 10^^. From these considerations we can already conclude that 
the proton decay rate in the current model can be bigger by a factor 100 compared to the 
the rate evaluated in four-dimensional SUSY SU(5) through dimension 6 operators. As 
we will see in the next section, we can obtain proton decay rates that are quite interesting 
for the next generation of experiments. Alternatively when considering option 0, where no 
flavour suppressions are present in the leading proton decay amplitudes, large values of Mc 
such as Mc 10^^ GeV allow to bring the proton lifetime above the experimental lower 
limit in eq. ((T)). Fig. 2 shows the predicted gauge coupling ay at the compactification 
scale Mc. As we can see from eq. ()36|) . this is an input to our computation and we evaluate 
it as the average among aj(Mc) [i = 1, 2, 3). We can see that au depends very mildly on 
the SUSY spectrum and on the treatment of heavy thresholds and is numerically close to 
1/25, the value in ordinary four- dimensional SUSY SU(5). The error bars in fig. 2 have 
the same meaning as in fig. 1. When computing proton decay amplitudes what matters 
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is the ratio au/M^, and we evaluate the errors on this quantity by fully accounting for 
correlations. 

Finally, fig. 3 displays the dependence of the ratio A/Mc on the SUSY spectrum, both 
for HN and CPRT cases. We can see that the average value of K/Mc is about 100, what 
needed for the dominance of the bulk symmetric SU(5) gauge coupling over the SU(5)- 
violating brane contribution in eq. ()48|1 . The ratio A/Mc is also directly related to the 
parameter e of eq. ()17|). which controls the hierarchies of Yukawa couplings. The range 
A/Mc ~ 100 is consistent, within the present approximations, with e ~ A^. 



5 Proton lifetime 



In order to calculate the proton lifetime, we have to translate the operators in ()28I36|37|) 
at quark level to those at hadron level using a perturbative chiral Lagrangian technique 
|33j . The aim is to evaluate the hadron matrix elements {PS\0\p), which describe the 
transition from the proton to a pseudoscalar meson via the three-quark operator O. The 
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Figure 4: Inverse decay rate of p ^ 7r°e+ versus the SUSY spectrum, in option 0. The solid 
line is the current experimental lower bound, from ref . P , and the dashed line represents 
the aimed- for future sensitivity, from ref. ^Hl- Error bars as in fig. CI 



various matrix elements are calculated from the basic element 



aPiUp 







(61) 



where Up denotes the proton spinor and a is evaluated by means of lattice QCD simulations. 
Here we will adopt 

|a| = 0.015 GeV^ . (62) 



''The statistical error of the lattice estimate is less than 10%, but the systematic error is much larger. 
For instance, within the lattice approach, ref. finds |q!| = 0.007 GeV'^. See also the discussion in ref. 
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In option 0, the dominant decay rate is given by: 

AlAljf^il + D + Ff + 



(63) 



Here mp, m^r denote the proton and pion masses respectively, and /tt = 130 MeV is 
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Figure 5: Inverse decay rate of p — >■ K versus the SUSY spectrum, in option I, for 
6 = (upper panel) and 6 = X'^ (lower panel). The solid line is the current experimental 
lower bound, from ref. pp. Error bars as in fig. ^ 

the pion decay constant; D and F are the symmetric and antisymmetric SU(3) reduced 
matrix elements for the axial-vector current and recent hyperon decays measurements j^ll 
give D = 0.80 and F = 0.46. The effect of lepton masses is neglected. The short- and 
long-distance renormalization factors As, encode the evolution from the GUT scale to 
the SUSY-breaking scale and the evolution from the SUSY-breaking scale to 1 GeV [SSj . 
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We use: 



Aj 



had) \ ' f a3{mc)\ f a^irub) _ ^ 43 



(64) 



with ^had = 1 GeV and 



A, 



' ai{mz) \ ^'""i / «2('"^z) \ ^''2 / d'iimz) \ ^''^ 



2.37 



(65) 



where the U(l) contribution is only approximate. In fig. E]we see the predicted inverse 
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Figure 6: Inverse decay rate of p ^ K^v (summed over all neutrino channels) versus the 
SUSY spectrum, in option I, for b = }? (upper panel) and 5 = (lower panel). The solid 
line is the current experimental lower bound, from ref . jl] , and the dashed line represents 
the aimed- for future sensitivity, from ref. ^U]- Error bars as in fig. ^ 

decay rate as a function of the SUSY spectrum. A general comment, that applies to all 
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the results concerning the rates in options 0, I and II, is that there is a huge theoretical 
uncertainty that spreads over many order of magnitudes. By far, the main source of this 
uncertainty is the compactification scale Mc, which is only know up to about two order 
of magnitudes, due to the unknown SU(5)-breaking brane contribution. Since the proton 
lifetime scales as M^, this corresponds to an uncertainty of more than eight order of 
magnitudes on the inverse rates. In this enormous range the probability is however almost 
uniform. From fig. ^ we see that most of the parameter space of the model have already 
been excluded by the existing experimental bound. Though the model is not entirely 
ruled out and the allowed region in parameter space will be almost completely explored 
by the next generation of experiments. The decay rates of option I and II for the different 
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(66) 
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Here mx denote the kaon mass and = 1.15 GeV is an average baryon mass according 
to contributions from diagrams with virtual S and A. In our numerical estimates the 
interference in the expression of the rate for K^v is constructive. The dominant proton 
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Figure 8: Inverse decay rate oi p K^v (summed over all neutrino channels) versus the 
SUSY spectrum, in option II, for 5 = \^ (upper panel) and 5 = (lower panel). The solid 
line is the current experimental lower bound, from ref . P , and the dashed line represents 
the aimed-for future sensitivity, from ref. ^Hl- Error bars as in fig. ^ 

decay rates are shown in fig. ()5I6|) and (pirUj). The results strongly depend on the location 
of Fi, either in the bulk (option I) or on the brane at ?/ = (option II). As we can see 
from Table 4, if Fi is a bulk field, then the dominant four-fermion operators are uus^ and 
udsu^^r and, consequently, the preferred proton decay channels are K^fi^ and K^u, whose 
rates are displayed in fig. ()5l6p . We see from fig. ()5I6|) that the possibilities of testing 
proton decay in option I are quite remote. If 6, which parametrizes our ignorance about 
the mixing angles, is equal to A^, the inverse rates are always larger than 10^^ yr, too long 
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to be observed in future planned facilities. When 6 = X^, only in some specific case the 
inverse rates can reach 10^^ yr, depending on the SUSY spectrum and on a favourable 
combination of the brane corrections 6^''^ . While the absolute rates are affected by a very 



p^7T*+v, Option II, 6-X' 



u 

>} 50 



O 

^ 45 
c 

-H 

g 40 

-H 
0) 

m 

-H 



35 



o 

±j 

o 30 
u 

CM 



la lb 23456789 SSSB 
SUSY spectra 



p^7T*+v, Option II, 6-X 



> 50 



o 

1^ 45 
C 

•H 

^ 40 

■H 
+J 

0) 

"1 35 



c 
o 

Jj 

o 30 

u 

CM 



la lb 2 3 4 5 6 7 
SUSY spectra 



9 SSSB 



Figure 9: Inverse decay rate of p tt^u (summed over all neutrino channels) versus the 
SUSY spectrum, in option II, for 6 = (upper panel) and 6 = X"^ (lower panel). The 
solid line is the current experimental lower bound, from ref. P]. Error bars as in fig. ^ 

large theoretical uncertainty, the predictions about the branching ratios are considerably 
more precise. The dependence on cancels out in the branching ratios and therefore the 
relative decay rates do not depend neither on the unknown SU(5)-breaking brane terms, 
nor on the SUSY spectrum. The uncertainty is dominated by the mixing matrices and is 
parametrized in our discussion by 6, which we let vary between and A^. In fig. [7|we 
see the BRs of the dominant channels in option I, K^u and K^fi'^, which are comparable 
within the estimated uncertainty. Substantially different is the option II, where Fi is a 
brane field like ^2^3. Here the preferred decay channel is K^u, as we can guess from Table 
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Figure 10: Inverse decay rate of p Tr'^e"'" (summed over all neutrino channels) versus the 
SUSY spectrum, in option II, for 6 = X'^ (upper panel) and 6 = X"^ (lower panel). The solid 
line is the current experimental lower bound, from ref . [T] , and the dashed line represents 
the aimed- for future sensitivity, from ref. fT^. Error bars as in fig. ^ 

6, showing that this final state has the smallest fiavour suppression. As we can see from 
fig. 8, a small portion of the parameter space of the model has already been excluded 
by the current experimental bound and, when 6 = X^, there will be good chances for the 
future facilities to discover proton decay in this channel, for almost all possible type of 
SUSY spectra considered in our analysis. Also if 6 = X^ and, consequently, the rates are 
more suppressed, there are several SUSY spectra which would allow a lifetime below 10^^ 
yr, at least for some combination of the brane kinetic parameters 6f'\ 

The rate for tt^u is down by approximately a factor A^, compared with the dominant 
one. Only for 6 = X^, for some specific SUSY spectra and in a favourable range of 5f , the 
predicted inverse rates are below 10^^ yr, which would however represent an improvement 
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over the current limit by more than two orders of magnitudes. 

The prediction for the decay channels containing a charged lepton are similar to those 
for TT'^u, but the expected future experimental sensitivity is quite better, especially in the 
^Og+ j-QQ(jg^ shown in fig. ^1 A modest portion of the parameter space, requiring 6 close 
to A'^, leads to inverse rates that are within the reach of future facilities. The branching 
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Figure 11: Branching ratios for option II. Dark- red (light-grey) histograms refer to 5 = A'^ 
(6 = X'). 

ratios for option II are displayed in fig. ^2 We can see the dominance of K~^u, followed 
by the tt^u and n^fi'^ channels. The other channels with a charged lepton, such as n^e^, 
K^e^ and K^fi'^ are comparable, but up to now future experimental prospects have been 
clearly worked out only for the vr^e"^ mode. 

6 Conclusion 

The decisive test of grand unification theories is proton decay. Minimal SUSY GUTs 
are strongly disfavoured today by the stringent experimental bounds on proton lifetime. 
Moreover such minimal realizations of the grand unification idea are also plagued by serious 
theoretical problems, first of all the doublet-triplet splitting problem, requiring a tuning by 
fourteen orders of magnitude. An elegant solution of the doublet-triplet splitting problem 
is offered by SUSY GUTs where the GUT symmetry is broken by the compactification 
of an extra dimension. There is a whole class of such GUTs where the leading baryon- 
violating operators arise from gauge vector boson exchange, have dimension six and scale 
as the inverse of the compactification mass Mc squared. 

We have evaluated form a next-to-leading analysis of gauge coupling unification 
including SUSY and GUT threshold corrections, two-loop running and SU(5)-breaking 
boundary terms. We have carefully estimated the uncertainties due to experimental errors 
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and poor theoretical knowledge, such as the ignorance about the SUSY breaking spectrum 
and about the details of the SU(5)-breaking brane dynamics. In view of the existing 
intrinsic theoretical uncertainties we have found a wide range of acceptable values for Mc, 
going from approximately 10^^ GeV to more than 10^^ GeV. 

The consequences on proton decay strongly depend on the features of the flavour 
sector, which can be considerably different from the standard, four-dimensional one. In 
the simplest case all matter fields have the same location in the extra space and proton 
decay proceeds mainly via the tt^c"*" channel, without any flavour suppression. Such a 
scenario is strongly disfavoured by now, since M^. is on average much smaller than the 
four-dimensional unification scale fa 2 x 10^^ GeV. However, due to the large uncertainties 
on Mc, it is not yet ruled out and future facilities can essentially complete the test of this 
model. 

When matter fields have no identical location in the extra dimension, fermion mass 
hierarchies and mixing angles are more naturally described. In our analysis we have 
considered two options, corresponding to "scmianarchy" (option I) and "anarchy" (option 
II) in the neutrino sector. In option I the suppression coming from fiavour mixing is 
too strong and proton decay is almost beyond the possibilities of the next generation of 
experiments. In option II the fiavour suppression of the channel can be compensated 
by the allowed smallness of resulting in a proton lifetime well within the reach of future 
facilities. 
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